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Accurate Monte Carlo critical exponents for Ising lattices 
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A careful Monte Carlo investigation of the phase transition very close to the critical point (T — > 
T c , H — > 0) in relatively large d = 3, s = 1/2 Ising lattices did produce critical exponents Pzd 
= 0.3126(4) ^ 5/16, 5^ = 0.1997(4) £ 1/5 and y iD = 1.253(4) £! 5/4. Our results indicate 
that, within experimental error, they are given by simple fractions corresponding to the linear 
interpolations between the respective two-dimensional (Onsager) and four-dimensional (mean field) 
critical exponents. An analysis of our inverse susceptibility data vs - l^ 1 — T c | shows that these 

data lead to a value of 7315 compatible with 7' = 7 and T c = 4.51152(12), while 7 values obtained 
recently by high and low temperature series expansions and renormalization group methods are not. 

PACS numbers: 64.60.-i, 64.60.Fr, 64.60.Cn. 



I. INTRODUCTION 

The Ising model [Q, rightly considered as the proto- 
type of statistical systems with non-classic power law 
critical behavior, has been extensively investigated for 
many years. Systems with short-range interactions dis- 
play Ising-like critical behavior, f.i. liquid-vapor, mul- 
ticomponent fluid mixtures, uniaxial magnets, etc, and 
there is a wealth of very accurate experimental informa- 
tion on these systems. 

To describe f.i. the behavior of a s = 1/2 uniaxial fer- 
romagnet near the critical point (T = T c , H = 0) two 
critical exponents, (3 for the spontaneous magnetization 
M S (T) and <5 _1 for the field dependence of the magneti- 
zation M C (H) at the critical temperature, determine ba- 
sically the critical behavior through M S (T) ~ \T C - Tf 
and M C {H) ~ H 1 / 5 . It is well known|| that for a two 
dimensional Ising lattice (d = 2) Onsager's solution gives 
fractional values for 02D = 1/8 and <5^ = 1/15. For a 
four dimensional Ising lattice (d = 4) on the other hand, 
the critical exponents are the mean field exponents j|, 
given also by fractional values, = 1/2 and = 1/3. 

It is a legitimate question to ask whether for a three- 
dimensional Ising lattice (d — 3), for which no general 
theoretical solution is available for the moment, the val- 
ues for and are rational fractions or not. In 
fact, almost forty years ago, Cyril Domb, one of the very 
pioneers in the then rapidly growing field of phase transi- 
tions, suggested that for three-dimensional Ising lattices 
the susceptibility critical exponent 7 = /3(<5 — 1) might 
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be given by the fractional value 73D = 5/4 = 1.25. Since 
then a tremendous amount of work (experimental, theo- 
retical and computational) has been performed with the 
aim to get ever more precise numerical characterizations 
of the phase transitions. Table I gives a representative 
sample j§ of numerical values §, |, @, |, |, [k], [ll], [l| for 
the exponents 73D and Psd obtained by various methods: 
high temperature expansion series, low temperature ex- 
pansion series, Monte Carlo simulations and field theoret- 
ical methods. The overall picture of the numerical values 
for 73 £> and (} 3 d is reasonably good, and they seem to 
favor 7 3 £> < 1.25 and fi^D > 0.3125, but, clearly, the un- 
certainties quoted in parentesis cannot be taken strictly 
at face value. 

In the present work we present results on critical expo- 
nents values based upon optimized accurate Monte Carlo 
calculations and we investigate to what extent the d = 3 
Ising exponents are compatible with the simple fractions 
interpolated between the fractional d — 2 Ising exponents 
and the, fractional too, d = 4 Ising exponents. In par- 
ticular we will use direct determinations of /?3D, from 
M s (T) data at H = 0, and of 5^, from M C (H) data at 
T = T c , as well as X _1 (^) vs T, including data both be- 
low T c (LT phase) and above T c (HT phase) , which allow 
us to make internal consistency checks, so that the actual 
value used for the critical temperature can be confirmed 
to be compatible with the scaling requirement 7' = 7 or 
not. 

In order to establish the reliability of the data and 
the propriety of the method of analysis used we will pro- 
ceed in two steps. First we will check data on large two- 
dimensional lattices, for which the fractional values of 
the exponents are known exactly, and then we will an- 
alyze data on large three-dimensional lattices for which 
the fractional values proposed are only educated guesses. 
Our data, therefore, can either lend support or leave un- 
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FIG. 1: 03d and 5^ (open symbols) obtained by interpo- 
lation between exact two-dimensional data from Onsager^] 
and four-dimensional data from mean fieldEJ theory (filled 
symbols). The crosses show current acceptedresultsfj]. In all 
cases, error bars are smaller than the symbol's size. Insets a) 
and b) are a blow-up close to the d = 3 region for /3 and (5 _1 , 
respectively. 

supported the fractional values proposed. Figure 1 shows 
the evolution of critical exponents with lattice dimension- 
ality d = 2, 3, 4,. . . for Ising systems. 



II. MONTE CARLO DETERMINATIONS OF 

M S (T) AND M C (T) 

Finite size scaling Monte Carlo simulations of phase 
transitions are known to be among the best and more 
effective techniques available to characterize the critical 
behavior of model systems such as Ising systems of any 
dimensionality. To determine critical exponents with the 
highest possible reliability it is desirable (i) to increase 
systematically the size L of the system till the results be- 
come practically independent of L within error bars; (ii) 
to make the interval (temperature or magnetic field) be- 
tween successive states in the vicinity of the critical point 
as small as possible, choosing wisely the full range (tem- 
perature of magnetic field) so that is not too large (time 
consuming) or too small (inconvenient of incomplete) to 
determine adequately the exponent in question; (iii) to 
insure that the time spent in the calculation at each point 
(temperature, magnetic field) is sufficient to arrive to 
true equilibrium. We did use two-dimensional lattices 
of 800 2 , 900 2 and 1000 2 spins, and three-dimensional lat- 
tices of 90 3 , 100 3 and 115 3 spins. The temperature in- 
tervals were of the order of 0.001, in units of T such that 
T c2D = 2.269185314213... and T c3D = 4.51152(12) for 
determinations of the spontaneous magnetization M s (T) , 
and field intervals of the order of 0.0005 for determina- 
tions of M C (H) at the critical isotherm. 

The number of Monte Carlo steps taken to insure equi- 
librium at each state was 50000 and the number of states 



considered in the partition function was 10000. No im- 
provement was detected by increasing the number of 
states. To get M S (T) a Wolff [[l4| algorithm was used, 
which is a generalization of the S wendsen- Wang [fl5| algo- 
rithm. To obtain M C (H), a somewhat novel but straight- 
forward, calculation not usually found in the literature, 
a standard Metropolis jl6| algorithm was used and it was 
proved to be quite adequate. 

Initially both periodic boundary conditions and free 
boundary conditions were used. After confirming that 
for very large lattices the difference between results ob- 
tained with either set of boundary conditions was quite 
negligible, we used subsequently periodic boundary con- 
ditions all the way. 

For the temperature scans we always began at T > T c 
and spent sufficient time at the beginning to make sure 
that thermal equilibrium was attained already far above 
the transition temperature. Nevertheless the statistical 
ups and downs in the residual M(T) data at T > T c , 
are much larger than the corresponding fluctuations at 
T < T c which become almost negligible and display a 
beautiful continuity all the way down in temperature. 
The log- log plots of M vs. \T-T C \ and M vs. H indicate 
clearly, first that we have been able to get really close to 
the critical point, more so than in other Monte Carlo cal- 
culations for which data in the literature are sufficiently 
explicit, and, second, that in the full range displayed for 
T or H corrections to scaling are invisible, which make 
the data especially apt to determine the exponents (3%d 
and Snn- 



III. MONTE CARLO RESULTS FOR 
TWO-DIMENSIONAL LATTICES 

Figure 2(a) gives a log-log plot of the spontaneous mag- 
netization as a function of temperature for a squared 
Ising lattice with 1000x1000 spins. Finite size effects 
show up as rounding at \T — T c \ — ► which occurs only 
at \T - T c \< 0.002 with T c = 2.269185314213... Correc- 
tions to scaling should appear at the other end of the 
temperature range examined, but are quite invisible in 
our data. The spontaneous magnetization data give di- 
rectly a value for f3 2D = 0.1242 ± 0.0008 1/8, as ex- 
pected. Figure 2(b) presents a log- log plot of the magne- 
tization as a function of field for the same squared Ising 
lattice with 1000x1000 spins at the critical isotherm. Fi- 
nite size effects begin to appear as incipient rounding at 
H < 0.0001 but are almost imperceptible. Corrections to 
scaling should appear at the other end of the field range 
investigated but are completely negligible in our data. 
These critical isotherm data, not previously investigated 
in depth, as far as we know, give directly = 0.06656 
± 0.00022 = 1/15, as expected. 

These results lend support to the expectation that 
carefully taken Monte Carlo data in large enough lat- 
tices taken at small enough temperature / field intervals 
are accurate enough to investigate whether critical expo- 
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TABLE I: Some estimated s — 1/2 critical exponents for d — 3 Ising simple cubic lattices 



Method Reference Year 7 /3 



High T expansions 




(2002) 
(1997) 


1.2368(10) 
1.2388(10) 


0.3243(30)" 
0.3278(13)" 


Low T expansions 




(1993) 
(1991) 


1.251(28) 
1.255(10) 


0.329(9) 
0.320(3) 


Monte Carlo 


X 


(1999) 

] (2000) 


1.2372(13)' 
1.255(18)' 


0.3269(5) 
0.325(5) 


Field Theoretical 




] (2001) 
] (1994) 


1.2403(8) 
1.258 


0.3257(5) 
0.336 



"Obtained using the scaling relation j3 = v(l + rf)/2 
'Obtained using the scaling relation 7 = (2 — rj)/v 



nents are given by simple fractions or not, at least in the 
case of two dimensions. 

A straightforward numerical analysis of the fractions 
compatible with the experimental results and the uncer- 
tainties quoted has been made. For d — 2 it can be seen 
that n/m with n < m compatible with the uncertainties, 
must go to m values very large, to > 256 for P2D and 
to > 360 for 5~f,- 



IV. MONTE CARLO RESULTS FOR 
THREE-DIMENSIONAL LATTICES 

Figure 3(a) gives a log- log plot of M S (T) vs. T for 
a simple cubic lattice with 115x115x115 spins. Lat- 
tices with 90x90x90 and 100x100x100 spins were in- 
vestigated and the results, of course slightly less accu- 
rate, where completely consistent with those obtained 
with the larger lattice. Again finite size effects appear at 
\T-T C \ — > 0, and they appear clearly at \T-T C \< 0.004 
but corrections to scaling are invisible in the temperature 
range investigated. For {3 3 d we get 

(i ZD = 0.3126(4) £ A (1) 
16 

Figure 3(b) produces a similar log- log plot of the mag- 
netization as a function of field for the same lattice with 
115x115x115 spins at T = T c = 4.51152(12). It is clear 
that both finite size effects and systematic departures 
from scaling are absent or imperceptible in the field range 
explored. These new, not previously directly investigated 
critical isotherm data, result in 

<y = 0.1997(4) - 1 (2) 

which differs somewhat from previous indirect es- 
timates of 63D summarized by A. Pelissetto and E. 
VicariQ as giving <5 3 £> = 4.789(2), equivalent to = 
0.2088(1). 



Using the numerical values for /? 3 £> and 5 3 p in Equa- 
tions (0) and (||) one gets 73D indirectly as 

J3D = PMho - 1) = 1.252(4) (3) 
by means of the scaling relation for 7 in terms of j3 and 

S. 

But one can get 7 directly from plots of susceptibility 
data at T > T c (j) as well as at T < T c (j') giving us 
the opportunity to check that the fundamental equality 
7 = 7', confirmed by renormalization group theory [Q, is 
fulfilled. 

Figure 3(c) gives log-log plots of direct Monte Carlo 
data on the susceptibility x~ 1 (T) as a function of \T — 
T c \ both above and below the critical temperature T c 
— 4.51152. The resulting directly determined value of 
gamma is 

73^=73^ = 1-253(4) (4) 

which is in excellent agreement with the forty-year-old 
prediction of Domb and Sykes ■ 

As mentioned before the following consistency check 
was made: changing smoothly the critical temperature 
value by increments (decrements) of 0.00002 we can ob- 
tain smooth changes in the effective values of the expo- 
nents 73D and "/' 3D to reproduce numerical values as low 
as 1.237 for 73c resulting in 7^ > 1.25, or as high as 
1.263 for 730, resulting in j' 3D < 1.25. With our Monte 
Carlo data, which in sufficiently wide ranges appear to 
be free of finite size effects as well as of corrections to 
scaling effects, only using the right T c — 4.51152, in per- 
fect agreement with the values quoted in the literature, 
the requirement 7 = 7' is duly fulfilled. 

Table II shows that according to our data for simple 
cubic s = 1/2 Ising lattices, simple fractions for (3 and (5 _1 
(the same procedure could be employed to check simple 
fractions for 7) compatible with the quoted uncertainties 
are given always by n/m fractions which are, either iden- 
tical to the interpolated values /3 3 _d = 5/16, d^f, = 1/5 
(l3D — 5/4) or much more complex fractions involving 
to > 128 (for (3) and to > 285 (for (S^ 1 ). 
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TABLE II: Monte Carlo s = 1/2 Ising critical exponents in a 
three-dimensional simple cubic lattice of 115 x 115 x 115 spins 



p 




n/m c 


0.31254±0.00029 


5/16, 


.., 40/128 (41/131) 






n/m c 


0.1997±0.0004 


1/5, . 


., 57/285 (57/286) 



fractional values compatible with the quoted uncertaintities up 
to rrimax = 128 for (3 and up to m max = 285 for 5 . More 
complex fractions (shown in parenthesis) are compatible with the 
uncertaintities only with denominators larger than 128 and 285, 
respectively. 



V. CONCLUDING REMARKS 

Our Monte Carlo data do not prove directly beyond 
doubt that three-dimensional Ising lattices are charac- 
terized by the fractional critical exponents @ 3D = 5/16 
and = 1/5 (and through the corresponding scaling 



relations, by the resulting critical exponents 7 = 5/4, 
a = 1/8, v = 5/8 and 77 = 0) but they support 
strongly and consistently this fractional values, in par- 
ticular the susceptibility data resulting in 73 £> = Y 3D = 
1.25 = 5/4. Our empirical data may be the basis for fu- 
ture well grounded theoretical arguments confirming the 
above fractional critical exponents. 
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FIG. 2: Two dimensional critical exponents (ho, S^d an d 720 
obtained by fits of Monte Carlo data of spontaneous magne- 
tization M S (H — 0) as a function of temperature T near 
T c = 2.269185314213, critical isotherm, M as a function of H 
at T = T c , and isothermal susceptibility x, respectively, in a 
system of 1000 x 1000 spins. Insets show the row data. 
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FIG. 3: Three-dimensional critical exponents 73D (a), fan 
(b) and 5^ (c) obtained by fits of Monte Carlo data of spon- 
taneous magnetization M s near T c = 4.51152 ± 0.00012, crit- 
ical isotherm and isothermal susceptibility Xj respectively, in 
a system of L 3 = 115x 115x 115 spins with periodic boundary 
conditions. Insets show the row data. 



